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We study triangular lattice spin-1/2 system with antiferromagnetic Heisenberg and ring exchanges 
using variational approach focusing on possible realization of spin liquid states. Trial spin liquid 
wave functions are obtained by Gutzwiller projection of fermionic mean field states and their en- 
ergetics is compared against magnetically ordered trial states. We find that in a range of the ring 
exchange coupling upon destroying the antiferromagnetic order, the best such spin liquid state is 
essentially a Gutzwiller-projected Fermi sea state. We propose this spin liquid with spinon Fermi 
surface as a candidate for the nonmagnetic insulating phase observed in the organic compound 
«:-(ET)2Cu2(CN)3, and describe some experimental consequences of this proposal. 



I. INTRODUCTION 

This paper reports a variational study of spin-1/2 
Heisenberg antiferromagnet with ring exchanges on a tri- 
angular lattice. One motivation for this study is the 
exact diagonalization work of LiMing et alX and Mis- 
guich et al^ on this system proposing that it realizes 
spin liquid states. We are particularly interested in spin 
liquids that may occur near the Heisenberg antiferromag- 
netic state. Multiple-electron exchanges are believed to 
be important near quantum melting and metal-insulator 
transitions. The specific model considered here may also 
be relevant for the description of a tentative spin liq- 
uid state observed in the quasi-two-dimensional organic 
compound k-(ET)2Cu2(CN)3^ which is close to metal- 
insulator transition. Imada et al^ studied appropriate 
Hubbard model on the triangular lattice and found an 
insulating regime with no spin order. The ring exchange 
spin model can be viewed as derived from the Hubbard 
model by a projective transformation, which is appropri- 
ate in the presence of the charge gap. 

The present work attempts to understand possible spin 
liquid states in the ring exchange model by examining 
candidate ground state wave functions. This is comple- 
mentary to the exact diagonalization studies, since know- 
ing the character of a candidate wave function can give 
significant intuition. 

The model Hamiltonian on the triangular lattice is, in 
the notation borrowed from Ref. 0, 
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The two-spin exchanges are between all nearest neigh- 
bors and reduce simply to Heisenberg interactions, P12 — 
P/2 — 25^1 • S'2 + 5. The four-spin "ring exchanges" are 
around all rhombi of the triangular lattice. 

In the following, we consider only antiferromagnetic 
coupling J2 > and positive J4 > 0; for brevity, we set 
J2 = 1. When J4 = 0, the system is the familiar Heisen- 
berg antiferromagnet on the triangular lattice and has 
a three-sublattice antiferromagnetic (AF) order. Exact 



AF 



spin liquid 

spin gap; many singlets 



-+ 
-0.1 



spin liquid 

no spin gap 



J4/J2 



FIG. 1: Phase diagram of th e model Eq. Q from exact diago- 
nalization study of Ref. I1I2I The magnetic order is destroyed 
for J4 > 0.07 — 0.1; a spin gap is observed in the regime 
0-1 ^ "^4 ^ 0.25, but also many singlets below the spin gap. 
The spin gap is decreasing for J4 > 0.175. 
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FIG. 2: Variational phase diagram for the Hamiltonian 
Eq. Q. The AF ordered variational state has the lowest 
energy for small J4, but becomes unstable for J4 > 0.14 
compared with the fermionic spin liquid states. One ex- 
ample of such spin liquid is the projected d^2_y2 + idxy su- 
perconductor ansatz, with the optimal variational parameter 
(A/t)var = 0.22,0.13,0.05,0.02 for J4 = 0.15,0.20,0.25,0.30 
respectively. Some other ansatze give very close optimal en- 
ergy, and the situation is particularly not clear near the AF 
state. But for J4 > 0.3, our best ansatze become essentially 
the projected Fermi sea state. We caution that for signifi- 
cantly larger J4 states with more complicated magnetic orders 
- e.g. with four-sublattice order - may enter the energetics 
competition,* which is not considered here. 



diagonalization study of Ref. proposes the phase dia- 
gram summarized in Fig. ^ The AF order is preserved 
for small J4 < 0.07 — 0.1, but is destroyed for larger 
J4 and a spin gap opens up. However, in the regime 
0.1 ^ Ja 0.25 reported in Ref . ilj, there are apparently 
many singlet excitations below the spin gap. Also, the 
spin gap starts to decrease for J4 > 0.175. 

In the exact diagonalization studies, it is hard to say 
which physical state is realized in the absence of clear sig- 
natures of some particular phase. The question of pos- 
sible spin liquid states is taken up here by considering 
variational spin liquid wave functions on the triangular 
lattice. Specifically, we consider one family of such states 
obtained by Gutzwiller-projecting singlet fermionic mean 
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field states?^ We determine the result of the competition 
with the AF ordered state by comparing against varia- 
tional wave functions with long range magnetic orderi 

The result of the variational study is summarized in 
Fig. El For small J4 < 0.14, the AF state is stable com- 
pared with the tried spin liquid states. For larger J4, 
we find spin liquid states that have lower variational en- 
ergy than the magnetically ordered state. For example, 
we find that projected superconductor ansatze perform 
well in the regime 0.14 ^ J4 < 0.3. More specifically, 
ansatze with anisotropic extended s-wave, dx2_y2, and 
dx2_y2 + idxy pairing patterns have very close optimal 
energies and much lower than the energy of the trial AF 
state. Unfortunately, we conclude that the present study 
is not sufficient to address the nature of the spin liquid in 
this regime, which we indicate with question marks in the 
figure. Our observation that the improvement in the trial 
energy is little sensitive to the specific pairing pattern 
may be an indication that the present restricted study 
cannot access the correct ground state in this regime. 

A more robust conclusion from our study of such spin 
liquids is that the best ansatze are close to the projected 
Fermi sea state and become more so for increasing J4. 
Thus, for J4 > 0.3 — 0.35 the variational A in our ansatze 
reduces to essentially zero (below few percent of the hop- 
ping amplitude), and the ground state is essentially the 
projected Fermi sea. 

The aptitude of the projected Fermi sea state can be 
intuitively understood as follows. We can view the ring 
exchange term with positive J4 as arising from the elec- 
tron hopping in the underlying Hubbard model (which 
we assume is in the insulating phase). Therefore, such 
ring exchange J4 > wants the fermions to be as "de- 
localized" as possible, and this "kinetic energy" is best 
satisfied in the simple hopping ansatz. A more formal 
mean field argument is given in Sec. IIIII 

We now discuss the indications of this study for the 
possible spin liquid state in k-(ET)2Cu2(CN)3. This ma- 
terial is close to the metal-insulator transition, so the 
role of the electron kinetic energy is clearly important. 
Based on the experience with the ring exchange model, 
we therefore propose that the projected Fermi sea state is 
a good candidate ground state close to the metallic phase. 
We verify this more explicitly by considering a model 
with ring exchanges obtained by a projective transfor- 
mation of the triangular lattice Hubbard model at order 
t^/U^. For the k-(ET)2Cu2(CN)3 compound we estimate 
J4/J2 ~ 0.3. The results are summarized in Fig.^] The 
work of Ref. on the triangular lattice Hubbard model 
can be interpreted as an elaborate numerical study build- 
ing up on free- fermion states, and hints some support to 
the proposed projected Fermi sea phase. 

The proposed picture has many physical consequences. 
We have a Fermi surface of spinous, and therefore expect 
no spin gap and finite spin susceptibility down to zero 
temperature consistent with the experimental observa- 
tions. An accurate treatment of the no-double-occupancy 
constraint and fiuctuations requires that the spinous are 



coupled to a fluctuating C/(l) gauge field. Such spinon- 
gauge field system has been studied extensively and is ex- 
pected to exhibit some unusual behavior 
For example, one expects a singular contribution to the 
specific heat Csing ^ T'^/'^ at low temperatures in two 
dimensions; the corresponding enhancement in "spin en- 
tropy" has concrete consequences for the phase bound- 
aries. 

The rest of the paper is organized as follows. In 
Sec.niwe specify the variational states considered in this 
work. In Sec. IIIII we seek qualitative understanding of 
the ring exchange energetics by considering a fermionic 
large- treatment of the ring exchange Hamiltonian. In 
Sec. IIVI we consider the connection with the triangular 
lattice Hubbard model and possible application to k- 
(ET)2Cu2(CN)3. In particular, we describe experimental 
signatures of the proposed spinon Fermi surface - gauge 
system. 



II. VARIATIONAL STATES AND ENERGETICS 

In this Section, we describe variational states used in 
the present work. Trial spin liquid states are constructed 
by Gutzwiller projection of fermionic mean field statesj^ 
We compare their energetics against AF ordered trial 
wave functions constructed using the approach of Huse 
and Elscr.^ 

Spin liquid trial states: The starting point here 
is the fermionic mean field treatment of the Heisenberg 
model. A recent and very detailed description can be 
found in Ref. 0- The setup for constructing trial wave 
functions is as follows. Each spin operator is written in 
terms of two fermions Cr-[ and c^j,, Sr = cl^Cr, with 
precisely one fermion per site. Heisenberg exchange in- 
teraction is written as a four-fermion interaction, which 
is then decoupled in the singlet channel. A convenient 
formulation of the mean field is to consider general spin 
rotation invariant trial Hamiltonian 
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with tr' r — ^rr' ' ^r'r — Aj.^. 



For each such trial Hamil- 
tonian we obtain the corresponding ground state. An 
SU{2) invariant formulation of the single occupancy con- 
straint is that the isospin operator Tr = ipl^^r is zero on 
each site; here ipri = c^i, i/'rX = c^,^, and r^'^''^ are Pauli 
matrices. In the mean field, we require that this con- 
straint is satisfied on average, which is achieved by tuning 
appropriate on-site terms. Going beyond the mean field, 
the physical spin wave function is obtained by projecting 
out double occupation of sites. 

Many such trial states can be constructed, but there is 
also a gauge redundancy in this construction. Here, one 
is helped considerably by the recently available classifica- 
tion scheme of X.-G. Wen^ that allows one to construct 
all possible such fermionic mean field states that lead to 
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physically distinct spin liquids with specified lattice sym- 
metries. 

We numerically evaluate the expectation values of the 
two-spin and four-spin exchanges in such states using 
standard determinantal wave function techniques (so- 
called variational Monte Carlo). We consider ansatze 
with different sets of lattice symmetries, with and with- 
out time reversal, but primarily we focus on the nearest- 
neighbor ansatze that respect upon projection the lattice 
translation symmetry. We then vary the parameters to 
optimize the trial energy. 



AF ordered trial states: We want to compare the 
ring exchange energetics of the spin liquid states with 
the energetics of the antiferromagnetically ordered states. 
For this purpose, we use the family of variational states 
considered by Huse and Elser^i which capture well the 
Heisenberg model energetics. Our primary goal here is to 
see how the AF state is disfavored by the ring exchanges. 
Since we are comparing with rather different states and 
are looking for the energy level crossing, we do not need 
to know the ground state energy very accurately, and the 
wave functions of Ref. should be sufficient to get rough 
idea of the ring exchange energetics in the AF state. For 
details on these wave functions and numerical evalua- 
tions, the reader is referred to the original paper. 



Variational results: We compared the trial energies 
of the AF ordered states and the fermionic spin liquid 
states, and the result is summarized in Fig. [3 For small 
J4, the ordered states have lower energy, but for J4 > 
0.14 the spin liquid states win. The optimal spin liquid 
ansatze have the following structure. The dominant part 
is the uniform triangular lattice hopping trr'-, and for 
J4 > 0.3 — 0.35 we essentially find the projected Fermi 
sea state. In the intermediate regime 0.14 < J4 ^ 0.3, we 
find that the trial energy is improved upon adding A^^' 
correlations into the mean field wave function. Somewhat 
perplexingly, we find that the result is not very sensitive 
to the specific "pairing" pattern. Thus, optimized wave 
functions with extended anisotropic s-wave, dx2_y2, and 
dx2_y2 -\- idxy pairing patterns have close energies. This 
may be an indication of an instability towards a state that 
cannot be captured in the context of the trial fermionic 
states. The situation is particularly inconclusive close 
to the AF phase, where several other trial states have 
competitive energies. 

In summary, we find that ring exchanges disfavor the 
AF ordered state compared with the fermionic spin liq- 
uid states, but our study is not conclusive as to which 
spin liquid state is realized when the transition happens. 
Away from the transition, we suggest that the optimal 
spin liquid state is the projected Fermi sea state. 



III. FERMIONIC LARGE N STUDY OF THE 
RING EXCHANGE ENERGETICS 

In this Section, we present a fermionic large N study 
of the ring exchange Hamiltonian. Here, natural "trial" 
states are pure hopping states, and this approach gives us 
some insight into their energetics. In particular, it shows 
how the ring exchanges favor the uniform hopping state, 
i.e. the projected Fermi sea state. The treatment below 
was suggested to the present author by T. Senthil. 

Consider the following generalization of the ring ex- 
change Hamiltonian Eq. (QJ to an SU (N) spin model 

HsU{N) = ^Xl("^l"^l'3)(4/3C2a) 
(12) 

+ J^Yl [(cLci/3)(4^C2^)(4^C35)(45C4a) + h.C. 
P 

We use conventional fermionic representation with N 
fermion flavors; spin states on each site are viewed as 
states of N/2 fermions, i.e. we have occupancy constraint 

cLcra = N/2 (3) 

for each site r. In the above, summation over repeated 
flavor indices is implied. Our generalization of the ex- 
change operators preserves the character of moving spins 
around a ring. For N — 2, this Hamiltonian reduces 
precisely to the spin- 1/2 Hamiltonian Eq. with 

J = 2J2, K = 8Ji. (4) 

A similar large N formulation was considered in a dif- 
ferent context in Ref. We also remark here that the 
general N formulation allows nontrivial exchanges involv- 
ing three spins. This is unlike the N = 2 case where such 
three-spin exchange reduces to a combination of two-spin 
exchanges. The three-spin exchanges can be easily in- 
cluded in the following analysis; to stay in line with the 
rest of the paper, we only consider the two-spin and four- 
spin exchanges. 

We formulate the large N procedure in the spirit of the 
variational approach. Consider a single-particle 'trial' 
Hamiltonian 

-Atrial ^ ^ (irr' ^ra^^' ^" h.C.^ ^ ^ f^r^ra^fOi- 

{rr') r 

We find the ground state and use it as a trial wave func- 
tion for the Hamiltonian Hgjj(^]^y In the mean field, the 
occupancy constraints are implemented on average by 
tuning the chemical potentials fir- The trial energy to 
leading order in 1 /N is given by 

^ = - \X12\^ - K^{X12X23X34X'11 + C.C.) , 

(12) P 

where xtr' = i'^l'^r') is the single-species expectation 
value. 



We now have to minimize -Emf over the possible trr' in 
the trial Hamiltonian. This leads to the following self- 
consistency conditions 



A ^trr' = JXrr' + ^ ^X23X34X41 > 

P=[1234] = [rr'34] 



(5) 



where the last sum is over all ring exchange plackets that 
contain the bond {rr') as one of the consecutive bonds. 
Also, we have explicitly indicated the fact that the trial 
energy does not depend on the absolute scale in the trial 
Hamiltonian but only on the relative pattern of trr' ■ 

We first make some general observations about this 
procedure. First of all, note that the self-consistency 
conditions imply that the optimal state can have non- 
zero trr' only on the bonds that have non-zero Jrr' or that 
appear in some ring exchange placket. For the triangular 
lattice model studied here, we then have to consider only 
nearest-neighbor trr' ■ Second, we see quite generally that 
the ring exchange contribution for a given placket has the 
form — ii'lxl* cos($p), where |xl is the geometric mean of 
the absolute values of Xrr' around the placket, while $p 
is the "flux" of the corresponding phase factors. Thus, 
the positive ring exchange wants to smear the fermions 
over the lattice with no fluxes. 

To be more precise, let us consider several simple trial 
states. Uniform flux state has flux (j) through each trian- 
gle. The expectation values Xrr' — {c\.iCr) have the same 
pattern of fluxes as the input i^r', and the trial energy 
per site is 



i?^ = -3J|x0p-6i^|X0|^cos(20) 



(6) 



since the flux through each rhombus is 2<j). Among such 
flux states, we flnd that for K < 2.76 J the best state has 
tt/2 flux through each placket (this state has the largest 
\x<p\)j while for K > 2.76J the best state has zero flux 
The numerical values of the energy per site in the two 
states can be obtained from 
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For large enough K/J the zero- flux state is stable against 
adding small flux (j) because |x0p/|xoP ~ 1 + 0.20^. 

We also considered so called dimer states such that 
non-zero trr' form nonoverlapping dimer covering of the 
lattice. These states break translational invariance, and 
any dimer covering produces such a state. It is well 
known that these states can have lower Heisenberg ex- 
change energy in the large N limit. This is because 
the occupied bonds attain the maximal expectation value 
I Xrr' I max and their contribution can be sufficient to pro- 
duce the lowest total energy. The energy per site in any 
dimer state is 

i^dimcr = -0.125J , (9) 

and is indeed the lowest energy for K = 0. However, the 
dimer states gain no ring exchange energy, and for K > 
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FIG. 3; Summary of the large study of the Hamiltonian 
Hsu(N)- a) Phase diagram from the mean field energy op- 
timization over translationally invariant states, b) Full opti- 
mization. 



9.9 J the zero flux state becomes the lowest energy state. 
Finally, the so called box states have identical two-spin 
exchange energy with the dimer states but also nontrivial 
fluxes and therefore do not enter the competition for the 
ground state for K > 0. 

We performed full optimization over trr' of the mean 
fleld energy considering possible unit cells with up to four 
sites, and found that the above simple states are indeed 
sufficient to describe the ground state in the large A'' 
limit: The optimal state is one of the dimer states for 
K < 9.9J and becomes the zero flux state for larger K. 
The complete study is summarized in Fig. |31 

To make better connection with the spin-1/2 system, 
we remark that we expect the dimer states to be energet- 
ically disfavored even for small K in the spin-1/2 case, 
e.g. compared with the flux states discussed above. This 
is because the Gutzwiller projection enhances local an- 
tiferromagnetic correlations more strongly in the trans- 
lationally invariant mean fleld states than in the dimer- 
ized statesii^ More quantitatively, the enhancement fac- 
tor for the Heisenberg energy is roughly (^yansi.mv. _ ^ 
for the translationally invariant states, while it is only 
^dimor _ 2 for the dimer states. Furthermore, we expect 
even stronger enhancements in the ring exchange energy 
upon the projection. Taking all this into account, we ex- 
pect the Fermi sea state to be favored for rather moderate 
J4/J2 in the spin-1/2 system. 

To conclude the mean field discussion, our main mes- 
sage is that the positive ring exchange dislikes the fluxes 
and wants to make the system as uniform as possible. 
This is best realized in the projected Fermi sea state. 

Going beyond the mean fleld, we obtain a theory of 
fermions coupled to a fluctuating gauge fleld (ao,a), 
where the temporal aq (r, r) enforce the local occupancy 
constraints while the spatial components represent the 
relevant fluctuations of trr'ir) « |t|e'"'-'-'(^). The corre- 
sponding continuum theory ( "relativistic electrodynam- 
ics in a metal") was studied in Refs. I Ml0llTlll]3l^l4lll5l 
and we will quote some results in the next section. 

In Appendix 1X1 we study long-distance properties of 
the Gutzwiller-projected wave function in some detail. 
As mentioned in the appendix, this wave function may 
be not sufficient to capture the long wavelength behav- 
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ior of the actual phase, since the projection treats only 
the Qq fluctuations, but does not include the fluctua- 
tions of Grr' , while the latter are crucial in the effective 
theoryiSiifliiiiiSiiiiiii^ This is pointing a possible limi- 
tation of the projected wave function approach for the 
spinon-gauge system. We still expect that the variational 
study of the previous section gets the crude energetics 
correctly in the ring exchange model. This is also what 
we expect from the mean field treatment, and leads us 
to propose the effective spinon-gauge theory. A finer nu- 
merical application likely requires more advanced tech- 
niques, perhaps in the spirit of Ref. for the triangular 
Hubbard model. It would be interesting for example to 
look for the 2kF signaturei^ in the more elaborate work 
of Ref. '4', which may be a more accurate realization of 
the spinon-gauge ground state. 

IV. APPLICATION TO POSSIBLE SPIN 
LIQUID STATE IN «:-(ET)2Cu2(CN)3 

We now discuss possible spin liquid state in the or- 
ganic compound k-(ET)2Cu2(CN)3 , which is insulating 
and shows no magnetic order down to the lowest exper- 
imental temperatures. It is believed^iliiS that the con- 
ducting layer of this material is well described by a single- 
band triangular lattice Hubbard model at half-filling with 
t/U ~ 1/8 and only small hopping anisotropy of about 
6%. 

Unlike the square lattice case, for the half-filled trian- 
gular lattice we expect a metallic phase for large enough 
t/U. Reference^ estimates the metal- insulator transition 
to occur at {t/U)Mi ^ 1/5, so the k-(ET)2Cu2(CN)3 ma- 
terial is on the insulating side. Using an elaborate nu- 
merical technique, Ref. |J finds a nonmagnetic insulator 
in this regime. We want to develop some picture of this 
state. 

The ideology we pursue here is that the insulating 
phase can be described by an effective spin model. Since 
the system is close to the metal-insulator transition, it 
is not enough to stop at two-spin exchange interactions. 
Starting with the Hubbard model, the effective Hamilto- 
nian to order /U^ was obtained in Ref. Specialized 
to the triangular lattice, the spin Hamiltonian reads 

iJeff = ^ring[J2, Ji] + ^ J" S, ' 5j + ^ j'" S, ■ (10) 

(m mm 

Here iJiing is the ring exchange Hamiltonian Eq. with 
J2 = (1 - 32tVC/2)2tVC/, J4 = 20tVt/3. The effec- 
tive Hamiltonian has additional Heisenberg exchanges 
J" = —l&t^/U^ between second neighbors (separated by 
a distance V3) and J'" = 4i'*/J7^ between third neigh- 
bors (separation 2 lattice spacings). Our grouping of the 
terms in the effective Hamiltonian is intended to make 
it look as close as possible to the ring exchange model 
studied in the previous sections. 

For the k-(ET)2Cu2(CN)3 compound, we estimate 
J4/J2 ~ 0.3, which puts the ring exchange model into 
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FIG. 4: Proposed phase diagram for the triangular lattice 
Hubbard model. The present study is based on the effective 
spin Hamiltonian Eq. llOt and applies only to the insulating 
regime expected for t/U < 1/5 from Ref. 0- Close to the 
metal-insulator transition, we propose the spin liquid state 
with spinon Fermi surface. For smaller t/U < 1/9, the best 
state is AF ordered. The k-(ET)2Cu2(CN)3 compound has 
t/U ~ 1/8 



the proposed spinon Fermi sea regime. Further neigh- 
bor interactions not included in the J2- model do not 
modify this result, even though J" and J'" are roughly 
of the same magnitude as J4. This stability is because 
the corresponding further neighbor spin correlations are 
small in the spin liquid regime. 

To proceed more systematically, we repeat the varia- 
tional study with the effective Hamiltonian H10|) . The 
resulting phase diagram is shown in Fig. 0] in terms of 
the Hubbard model parameter t/U . From this study, we 
propose that the insulating ground state is the antiferro- 
magnet for t/U < 1/9 (this corresponds roughly to the 
ring exchange parameter J4/ J2 ~ 0.2 — 0.25). For larger 
t/U , our best trial state is essentially the projected Fermi 
sea state, and the variational A (which can be used to 
improve the trial energy slightly) is small already at the 
transition from the AF state. In the same figure, we also 
indicate the metallic phase expected for t/U >^ 1/5. 

It should be emphasized that we do not treat either 
Hamiltonian Eq. H10|) or (^) as more realistic or less realis- 
tic, particularly since we are dealing with the system near 
the metal-insulator transition. The above variational 
study with Hcs is presented primarily to illustrate that 
our results are not destabilized by making the Hamilto- 
nian 'more realistic'. We expect that our main prediction 
for the spin liquid state close to the metal-insulator tran- 
sition is robust, since the proposed Gutzwiller-projected 
Fermi sea state is even more favored by including further 
effects of the electron kinetic energy. Also, the results 
of Ref. i give us some indication on the stability of the 
proposed state, since that study is building up on free- 
fermion states. 



A. Physical properties in the spin liquid phase 
with spinon Fermi surface 

The effective description of the proposed phase has 
spinon Fermi sea coupled to a dynamically generated 
gauge field. It has been argued^i^^i^^'^^i^^i^'^i^^ that this 
spinon-gauge system is described by a nontrivial fixed 
point and shows unusual behavior, which can be tested in 
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experiments. Below, we list some thermodynamic prop- 
erties of this Mott insulator. This phase is in some sense 
the closest one can get to the Fermi liquid while remain- 
ing a charge insulator, and shares some properties with 
the metal due to the presence of the spinon Fermi surface, 
but also has some 'non-Fermi liquid' properties. 

Thus, spin susceptibility is expected to approach a 
constant as temperature T goes to zero: 

Xspin(r 0) /i|i/o . (11) 

This is a consequence of having gapless spinon excitations 
over the entire Fermi surface and is in fact observed in k- 
(ET)2Cu2(CN)3 ^ Here, vq is the density of states at the 
"Fermi surface" in the spinon band structure determined 
by the spinon "hopping amplitude" tspinon- The latter is 
set by the Heisenberg exchange energy igpinon ~ J and 
is different from the bare electron hopping amplitude td 
(remember that J ~ t^^/U). For the triangular lattice at 
half-filling, we have = 0.28/isninon per triangular lat- 
tice site and including spin. Ref.^ reports x = 2.9 x lO^'* 
emu/mol at low temperatures, from which we estimate 
^spinon ~ 350K. This compares favorably with the re- 
ported magnitude of the Heisenberg exchange J = 250K. 

Specific heat, on the other hand, is expected to show 
non-Fermi liquid behavior 

C-kBV,tlll^^{kBTf/^. (12) 

This is written to contrast with the Fermi liquid ~ T 
behavior, and means that the spin entropy in this charge 
insulator is in fact larger than in the metallic state at 
low temperature. This is very different from the anti- 
ferromagnet or gapped spin liquid insulators which have 
low spin entropy. In particular, the finite temperature 
first-order transition line between the proposed spin liq- 
uid and the metallic state is expected to bend towards 
the metallic state with increasing temperaturei^ 

mi(T)-PMi(o)~r'/^ (13) 

In the last formula, p is an applied pressure which drives 
the insulator to metal transition. ^'^^ This tendency is ac- 
tually observed in the k-(ET)2Cu2(CN)3 materiali^i 



V. CONCLUSIONS 

In summary, we considered the spin-1/2 ring exchange 
model on the triangular lattice from the variational per- 
spective and identified the instability of the antiferro- 
magnetically ordered state towards spin liquid state in 
the regime of moderate ring exchange couplings. Our 
best trial states become the Gutzwiller-projected Fermi 
sea state for larger J4. Despite the limitations of the vari- 
ational approach, it is hoped that the present work may 
give complimentary information and useful guidance for 
understanding the exact diagonalization results. 



We also studied the effective spin Hamiltonian appro- 
priate for describing charge insulator states of the trian- 
gular lattice Hubbard model. The effective Hamiltonian 
includes Heisenberg exchanges as well as ring exchanges, 
and so is close to the considered ring exchange model. 
The study is motivated by the tentative spin liquid state 
in the k-(ET)2Cu2(CN)3 compound, which is modeled 
by the triangular lattice Hubbard model in the vicin- 
ity of the metal-insulator transition. We find that upon 
including the ring exchanges but well in the insulating 
regime, the antiferromagnet gives way to the spin liquid 
state which is essentially the projected Fermi sea state. In 
view of this finding, we propose that the effective descrip- 
tion of the nonmagnetic insulator phase has Fermi sea of 
spinons coupled to the dynamically generated gauge field. 
This spin liquid phase features a number of unusual prop- 
erties which can be looked for in experiments. It would 
be very exciting if this remarkable state is indeed realized 
in the k-(ET)2Cu2(CN)3 material. 
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APPENDIX A: PROPERTIES OF THE 
PROJECTED FERMI SEA WAVE FUNCTION 

We describe some properties of the projected Fermi sea 
state. Figure 13 shows spin correlations in the projected 
wave function and also in the free fermion state before 
the projection. In the free fermion state, the spin correla- 
tion behaves as — cos^(A:Fr — 37r/4)/r^ at large distances, 
which oscillates with the wave vector 2kp while always 
staying negative. To facilitate the comparison. Fig. |S1 
shows the mean field correlations in the specific finite sys- 
tem (the finite size effects are fairly large because of the 
gaplessness over the Fermi surface). We observe that the 
effect of the projection is not strong: For the range stud- 
ied, the mean field result multiplied by the Gutzwiller 
enhancement factor (?j = 4 gives a reasonable approxi- 
mation for the actual correlationiiS After the projection, 
the correlation function now swings to positive values as 
well, but the overall magnitude is roughly captured by 
the simple renormalization factor. 

We also studied spin chirality correlations (not shown) , 
and found that these are very small beyond few lattice 
spacings. 

The effective theory of the proposed phase has Fermi 
sea of spinons coupled to a dynamically generated gauge 
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FIG. 5: Spin correlation in the projected Fermi sea state. 
Measurements are done on a 24 x 24 triangular lattice. The 
mean field wave function is constructed for periodic boundary 
conditions and excluding the zero momentum single particle 
state in order to avoid Fermi surface points while satisfying 
the lattice rotation symmetry for the finite system (this does 
not affect the long distance properties of the wave function 
which is our focus here) . Note the oscillating character of the 
correlation (with the period ~ 2n/{2kF), kp ~ 2.69). Also 
note that the renormalized mean field roughly reproduces the 
overall magnitude of the correlations. 



fieldi2ii2iiLi2iiii4 The measured spin correlations in the 
projected wave function represent some puzzle in this re- 
spect: Ref. 13 predicts that the spin structure factor is 
singularly enhanced near 2k p in the spinon-gauge sys- 
tem. We find that in the projected wave function the 
structure factor remains finite throughout the Brillouin 
zone and that the overall rate of decay of spin correlations 
is roughly the same as in the free fermion state. One pos- 
sible source of this difference is that the projected wave 
function has fixed trr' and therefore does not include the 
fluctuations of the spatial components of the gauge field; 
only the temporal component is "included" by the projec- 
tion. This is a limitation of the projected wave function 
approach for the spinon-gauge system. 
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